Lagrangian Particle Statistics in Turbulent Flows from a Simple Vortex Model 
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The statistics of Lagrangian particles in turbulent flows is considered in the framework of a simple 
vortex model. Here, the turbulent velocity field is represented by a temporal sequence of Burgers 
vortices of different circulation, strain, and orientation. Based on suitable assumptions about the 
vortices' statistical properties, the statistics of the velocity increments is derived. In particular, the 
origin and nature of small-scale intermittency in this model is investigated both numerically and 
analytically. 
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I. INTRODUCTION 

Our understanding of the spatio-temporal properties 
of turbulent flows is still fragmentary [U, 0, Q . In recent 
years, however, significant progress could be achieved by 
focusing on the dynamics of Lagrangian particles (see, 
e.g., Ref. Q and references therein). Especially labo- 
ratory experiments have given useful insight to impor- 
tant Lagrangian statistical properties. For example, the 
statistics of Lagrangian velocity increments 0,1^ and 
those of the turbulent acceleration of a fluid particle 0, Q 
have been measured. Moreover, direct numerical simu- 
lations on large-scale computers have been performed [Sj 
which allow to study particle trajectories and their sta- 
tistical properties from yet another point of view. While 
all results obtained from laboratory and computer ex- 
periments seem to be consistent so far, the underlying 
physical mechanisms are still subject to discussion. In 
particular, it has not yet been possible to arrive at a theo- 
retical derivation of the single particle velocity increment 
distribution. 

There are hints, though, that certain turbulent struc- 
tures are largely responsible for the statistical properties 
of Lagrangian particles. It has been shown by direct nu- 
merical simulations |p, [lo| and by thorough experimen- 
tal investigations that the turbulent field contains 
elongated vortex filaments which can be interpreted as 
Burgers vortices. The latter are well-known solutions of 
the Navier-Stokes equation [l^. In fact, these vortices 
have been called the sinews of turbulence by Moffat and 
coworkers 13|. Hatakeyama and Kambe recently mod- 



elled turbulent fields by a random arrangement of such 
Burgers vortices Interestingly, through this proce- 

dure they were able to reproduce the multifractal scaling 
behavior of the longitudinal structure functions Snir). 

In the present paper, we shall employ a similar ap- 
proach. More precisely, we will consider the Lagrangian 



statistics of a tracer particle in a turbulent velocity field 
which is represented by a temporal sequence of Burg- 
ers vortices of different circulation, strain, and orienta- 
tion. Such a model is motivated by the fact that Navier- 
Stokes turbulence is known to create strong vortex fila- 
ments which tend to dominate the time evolution of La- 
grangian particles. This view is supported, e.g., by the 
recent work of Biferale and coworkers [1, [3] who detect 
so-called "vortex trapping events" of tracer particles in 
their direct numerical simulations. Hence we are led to 
consider the path of a Lagrangian particle in the field of 
a single Burgers vortex. After a certain lifetime T/, this 
vortex decays and is replaced by another vortex which 
differs in circulation F, strain-parameter a, and orienta- 
tion. Then, the process starts anew. Making suitable 
assumptions about the vortices' statistical properties, we 
will be able to determine the statistics of the velocity 
increments in the framework of our model this way. In 
particular, we will be able to investigate the origin and 
nature of small-scale intermittency both numerically and 
analytically. 

The remainder of the present article is structured as 
follows. Starting from the exact solution of a particle 
trajectory in a Burgers vortex, we will investigate the 
functional structure of the probability density function 
(pdf) of the velocity increment. We then will present 
results of a Lagrangian particle evolving through a tem- 
poral sequence of trapping events. There especially the 
evolution of the pdf 's of the velocity increments is con- 
sidered. 



II. LAGRANGIAN PARTICLE IN A SINGLE 
BURGERS VORTEX 

The velocity field of a single Burgers vortex with cir- 
culation F in a strain field Ust(x, t) = [—fa;, — fy,az] is 
given by 
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u{x, t) = u,t{x, t)+e,x Br — [I - e-"'^'/(4"^)] (1) 
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where v denotes the kinematic viscosity of the fluid, F de- 
notes the circulation of the vortex and r = \J + y^- e-z 
and denote the unit vectors in z- and radial direction, 
respectively. Switching now to the Lagrangian frame, the 
evolution of a Lagrangian particle starting from an initial 
position cco is determined by 



dX 



dt 
<PX_ 
dt^ 



(2) 



{xo,t) = [-Vpix,t) + iyAuix,t)]^^x^^^^yi3) 



Due to the axial symmetry of the problem, we seek for 
solutions of the form 



x{t)=r{t) eriifiit)) + z{t) e. 



(4) 



where r{t), cp{t) and z{t) have to obey the following set 
of differential equations: 



r = r 
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with zi„ = - e-'"-'/(4!^)l. The solutions for two of 
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these components obviously read 



z{t) 
r{t) 
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So the influence of the z-dynamics increases exponen- 
tially. Since we want to focus on the importance of the 
oscillatory motion around the axis of the vortex, we have 
to study the case of low straining, i.e. the factor at has to 
be small compared to the circulation F. Plugging the so- 
lution for r{t) into the ODE for the azimuthal component 
one obtains 



iP{t) = 



2^r2 



1 -e" 



(10) 



The solution of this ODE is thereby reduced to an inte- 
gration. However this integral cannot be solved explicitly, 
and hence we want to restrict our following calculations 
to the limiting cases where the particle is far away from 
the vortex's axis or near the core. For the farfield case the 
point-vortex approximation holds, leading to an angular 
velocity of 



(11) 



The angular velocity of the particle decays like lead- 

''o 

ing to a differential rotation. In case of low straining 
the exponential can be approximated by unity. In this 
approximation the solution is 



(12) 



Note that due to this approximation the angular velocity 
becomes independent of the strain-parameter a. 
The second limiting case is where the particle is near the 
viscous core of the vortex. In this case we can expand the 
exponential into a series. The first order approximation 
then is 



which leads to 
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Lp{t) = - — i = uj{T,a)t. 



(13) 
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This means that near the viscous core of the vortex the 
motion is given by a rigid body rotation determined by 
the parameters a, F and v. In order to clarify the func- 
tional structure of the following results for the velocity 
increment distribution we will first calculate the simple 
case, where only one vortex is involved. 
Statistical observations on turbulence often focus on the 
velocityincrements Vx{t) with a time lag r (see for ex- 
ample 0) and on the turbulent acceleration [l^. For the 
case of small t the evolution of a particle in a turbulent 
flow is dominated by the nearest vortex filament. So the 
statistics for small r can be derived directly from the dy- 
namical equations of a single vortex. 
Hence we have to write down the discussed solution in 
cartesian coordinates. The x-component of the position 
of the Lagrangian particle is given by 



x{t) — roe 2* cos fit). 



(15) 



Intermittency in the velocity signal is often said to be 
caused by the strong accelerations of a vortex filament 
■ These accelerations originate from the oscillation of 
the particle round the axis of the vortex [18j . We now 
neglect the radial transport of the particle in order to 
focus on these rotational accelerations. This approxima- 
tion holds whenever the product at is small compared 
to the vorticity determined by F. As we want to focus 
on the impact of strong vortex filaments on the statis- 
tics this is a reasonable approximation. Neglecting the 
straining the x-component of the velocity reads 
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(16) 



with uj being the oscillation frequency given by the two 
limiting cases discussed above. This leads to a simple 
equation for the velocity increments, 

Vx{t) = Uxit + r) ~ Ux{t) 

— —rouj {sin u!(t + t) ~ sin ut) . (17) 

Some simple trigonometric relations turn this expression 
to 
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(18) 
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with a phase ^p which depends on the initial conditions. 
Note that the velocity increments oscillate like the ve- 
locity components themselves, being modulated by an 
amplitude which depends on the parameters of the vor- 
tex as well as on the time lag r. The next aim is to 
deduce the corresponding probability density function. 
Instead of taking a time average, we average over the 
phase a = cot + ip, which can be assumed to be uniformly 
distributed, 

1 f^- 

f{va;) = ^ / daS{vx ~ Asma) 

1 f^'^ 1 Vx 

da — r 5{a — arcsin —^). (19) 
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Whenever < A the solution reads 
1 



1 



(20) 



The dependence on the parameters F and a is absorbed 
into the frequency to. Noting that the violation oivx < A 
leads to imaginary values of the probability, the restric- 
tion becomes redundant when expressing the solution as 
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Figure 1: Trajectory of a Lagrangian particle, t is given in 
multiples of t^. 



5ft 



(21) 



In order to obtain our result we have to average over 
the vorticity parameter F, the strain-parameter a and 
the initial position of the particle, vq respectively. Let 
/(F), g(a) and h{ro) be the corresponding distributions 
of the parameters. (This notation implies statistical in- 
dependence, a restriction, which has not to be made). 
The probability distribution of the velocity increments is 
then given by 



/(F)ff(a)/t(ro) dFdarpdr o 



(22) 



This expression explicitly reveals the dependence of the 
probability density functions on the distribution of the 
vorticity parameter F. Thereby a physical connection 
between the velocity increment distribution and the 
physical parameters of the flow is established. 
The structure of the functional form of the pdf reveals 
that the shape of the pdf is given by a superposition 



of functions of the form 



weighted by factors 



from the physical parameter distributions. The complex 
interplay of the parameter distribution then leads to the 
fat-tail pdf's, as observed in the following. 



III. LAGRANGIAN PARTICLE IN A 
SEQUENCE OF VORTICES 

As discussed above, we model the evolution of a La- 
grangian particle in a turbulent flow by a temporal se- 
quence of randomly oriented coherent structures. In par- 
ticular a particle is exposed to the velocity field of a sin- 
gle vortex for a given lifetime T/. In our model T; is 
given by the typical lifetime of an eddy in the dissipa- 
tive range. According to this lifetime is suggested 
to be of the order IOt^ (r^ — {e/v)^ denotes the Kol- 
mogorov timescale). In real turbulent flows viscosity or 
vortex merging processes then force this vortex to de- 
cay. This is crudely modelled by simply turning the vor- 
tex off. Subsequently another vortex is switched on at 
a fixed distance. This distance is given by the typical 
spatial density of strong vortex filaments in a turbulent 
flow (see for example [Sj], chapter 7.3). The orientation of 
the vortex is chosen randomly. As a result the statistics 
become isotropic. 

In order to render the model more physical, some prop- 
erties of turbulent flows have to be taken into account. 
Special interest has to be put on the statistical prop- 
erties of the strain-parameter a and the circulation F. 
A hint on the distribution of F is given by the authors 
of 0, who measured the distribution of the vortices' 
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strengths directly from numerical simulations. In our 
model this is taken into account by modelling the dis- 
tribution of r with a log-normal distribution. Also the 
strain-parameter a is chosen randomly. This means, that 
in our model vortices of different radii appear, as the 
strain-parameter a directly affects the radius of the Burg- 
ers vortex. However, this distribution cannot be equally 
distributed, because this would cause vortices of any pos- 
sible size to appear. Therefore a peaked distribution is 
needed. The expectation value of this distribution is re- 
lated to the typical size of the smallest vortex filaments. 
We choose this distribution to be Gaussian with standard 
deviation ct^ = 1- This choice is somewhat arbitrary, but 
it was ensured that the numerical results vary only little 
over a wide range of possible values of aa- The expec- 
tation value of the distribution can be derived from the 
definition of the radius of the Burgers vortex. This typ- 
ical scale is defined as = (^)^. If we demand rs 
to be of the order of IO77, a value which conforms to 

the findings of the expectation value of the strain- 

1 

parameter a is given asa= ^{f:)^ = 557"' '^^ere e 
denotes the energy dissipation of the system. Note that 
by this relation a connection between typical scales of the 
velocity field of the Burgers vortex and physical proper- 
ties of the flow is established. Moreover, the straining 
thereby is related to the Kolmogorov time scale r,,. A 
typical value of the circulation has to be determined in 
the same manner. We start with the well-known relation 
for the Reynolds number on the dissipative scale of the 
flow, Remiss = Ri 1. The typical velocity can be 
approximated by the velocity field of the Burgers vortex 
on the scale 77, 
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MV) = ^ ( 1 



Far; 
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T^- (23) 
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This yields F = — {ei^)^ = 2007rz^, which determines a 
typical value for F. The distribution of F is assumed to 
be log-normal according to the observations of [9|]. 
Now all properties of the model are specified and we are 
able to perform a numerical simulation. 
Figure [T] shows a typical particle path. It reveals that the 
path is composed of a sequence of vortex trapping events, 
each of them characterized by a different circulation F. 
Events with a strong circulation can easily be identified 
by their rapid oscillation of the a;-component, whereas 
events with small values of F do not cause the particle to 
perform a complete oscillation. The switching between 
vortices renders the velocities discontinuously. This dis- 
continuity can be eliminated using a low pass filter. A 
lowpass filtering of the velocity signal with a Butterworth 
filter was performed and the resulting statistics was com- 
pared to the unfiltered one. It was ensured that this effect 
has no significant impact on the statistics. 
One of the central results of our model is depicted in fig- 
ure [2] which shows the probability distribution of the ve- 
locity increments. These exhibit a similar behavior to the 
velocity increment distributions observed in experiments 




Figure 2: Logarithmic plot of the pdf's of the velocity in- 
crements. Time lags Te{ (0.98, 1.97, 3.94, 7.87, 15. 74)r^}. The 
pdf's are normalized to cr = 1 and shifted vertically for view- 
ing convenience. 




Figure 3: Kurtosis K{t) = 
let shows a log-plot of K{t 
exponential decay. 



The in- 



. The kurtosis shows an almost 



or direct numerical simulation. For small time lags t the 
pdf's show fat tails, whereas a nearly Gaussian shape 
arises for large time lags r. The physical explanation is 
straight forward; while the functional form for the short 
time lags is almost solely determined by a single vortex, 
the transition to a more Gaussian functional form results 
from an increasing statistical independence as more vor- 
tices are involved. The transition from one functional 
form to the other is characterized in more detail by the 
kurtosis, which is shown in figure [31 This figure clearly 
indicates the transition from a highly intermittent dis- 
tribution to a more Gaussian one. Since the subsequent 
trapping events are statistically independent, the velocity 
signal has to decorrelate as a function of the time delay 
r. We define the velocity autocorrelation function as 



C{t) 



(24) 



Figure 2] confirms this view; the correlation decreases 
almost exponentially. The slight anticorrelations might 
tribute to the events with strong circulations, where the 
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Figure 4: Autocorrelation function of the velocity u^;, r is 
given in multiples of t,,. The inlet shows a log- plot of the 
vertically shifted autocorrelation function. 

particle is able to perform one or more complete oscil- 
lation. The velocity signal decorrelates after a time of 
T = lOr,,, which is the lifetime of our vortices. A second 
effect which cannot be clearly separated is dephasing. As 
the particle moves inward it constantly alters its angular 
velocity. It can be shown easily that this also causes the 
autocorrelation function to decrease. 



IV. SUMMARY 

To conclude, we modelled the particle evolution in a 
turbulent flow as a temporal sequence of Burgers vor- 
tices. We are able to calculate the path of a Lagrangian 
particle in an isolated vortex. From this point we derived 
the functional form of the velocity increment pdf. This 
calculation reveals the dependence of the structure of the 



velocity increment distribution on the physical parame- 
ters of the vortex filament like strain, circulation etc. 
A numerical implementation of the model was used to 
investigate some typical Lagrangian properties of turbu- 
lent flows. As a central Lagrangian observable we fo- 
cused on the velocity increment distribution finding that 
our model qualitatively resembles the velocity increment 
statistics in real turbulent flows. A further investigation 
of the corresponding kurtosis clearly indicated intermit- 
tent characteristics. 

Additionally the velocity autocorrelation function was in- 
vestigated revealing an almost exponential decrease and 
even slight anticorrelations. Consistent with the frame- 
work of the model is the fact that the autocorrelation 
function vanishes for times longer than the lifetime of a 
single vortex. 

The physical interpretation of the observed results is 
quite straightforward. A Lagrangian particle encounter- 
ing subsequent trapping events shows intermittent veloc- 
ity increment distributions accountable to a proper sta- 
tistical ensemble of Burgers vortices. For small time lags 
the fat-tail statistics originate from the superposition of 
the dynamics in a single vortex. The transition to a more 
Gaussian behavior is caused by an increasing mixing of 
two subsequent vortex trapping events which are statisti- 
cally independent. This observation is supported by the 
functional form of the velocity autocorrelation function. 
This model qualitatively reproduces some of the typical 
Lagrangian statistics. It is important to note that this is 
achieved by a temporal composition of exact solutions of 
the Navier-Stokes equation. This is in contrast to many 
other models, which apply stochastic equations for the 
particle evolution. Due to the simplicity of the model 
our results can be interpreted in physical terms and shed 
light on the connection between dynamical aspects in tur- 
bulent flows and corresponding statistical properties. 
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